From Quantum to Classical: the Quantum State Diffusion model 
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1 Introduction 



The vagueness of the quantum-classical border in physics makes some physicists uncertain about 
| the foundation of their science. Others, by contrast, feel that the accuracy of the verified 
predictions makes these foundations especially solid. This longtime state of affairs is changing 
in an interesting way. Recent progress in mesoscopic physics, and the so-called "new science 
of nanotechnology," open new domains to be explored by theoreticians and experimentalists 
alike. In this contribution we will exploit numerical and conceptual tools that were developed 
to explore this new domain of mesoscopic physics, and apply them to systems that are close to 
^ ', the quantum-classical border. 

Central to classical mechanics are the concepts of trajectories and the phase space structures 
associated with them. Hence, we would like to illustrate our results with figures showing these 
trajectories and structures as they appear in our models. These figures, and the fact that they 
can be computed without any supplementary parameters, constitute the main result of this 
contribution. 

As examples we could use several kinds of oscillators. The first, and the simplest, is the 
harmonic oscillator with its regular orbits; but this is boring! Next, systems that exhibit classical 
chaos, such as the Kicked Anharmonic OScillator (KAOS) ||, |3], [| and the double-well Duffing 
oscillator || §]. In both cases the structure of the strange attractors emerges out of quantum 
cloudiness when the quantum-classical border is approached. Another example worth mentioning 
is a nonlinear oscillator exhibiting hysteresis |?J , the hysteresis curve appearing when the system 
becomes more classical. 

The feature which characterizes all these systems is that they are open; they interact signifi- 
cantly with some kind of external environment. One cannot in general describe an open system 
by a single state vector; one requires a density matrix, essentially giving classical probabilities 
for the system to be in given quantum states. In the limit of Markovian systems, this density 
matrix evolves according to a linear master equation, as we shall see below. 
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For many reasons, we would like to get away from the density operator description, and 
these new tools allow us to do so. These tools require no additional parameters, only a simple 
idea: describe open quantum systems by state vectors, as for closed systems. Unlike closed 
systems, their evolution is not unitary — in fact, it is not linear and not deterministic (two 
general characteristics of open systems). The evolution is chosen such that when the pure state 
defined by the (normalized) state vector is averaged over the randomness, the usual density 
matrix description is recovered at all times. The numerical advantage of this approach stems 
from replacing an n x n matrix by an n component vector. The conceptual advantage lies in 
the possibility to have individual states in the theory, as in actual experiments, and to compute 
mean values as averages of individual outcomes, again as in actual experiments. 

This approach shares with the decoherence approach the idea that the environment plays a 
crucial role in breaking the unitary (boring) evolution and reducing the state § [|. However, 
here the evolution provides classical probabilities for quantum amplitudes. Our approach also 
shares some similarities with the GRW spontaneous localization model J10H , in that the evolution 
equations are stochastic, but preserve pure states. But here there are no ad hoc parameters. 

2 Quantum State Diffusion and Quantum Jumps 

Markovian open quantum systems are usually described by a master equation: 

P = ~i[H,p] +J2( L "ipLln- \{ L ln L m,p}\ (1) 

where p is the density matrix for the system, H its Hamiltonian, and the linear operators L m 
describe the effects of the environment. In an unraveling of such master equations, one describes 
the system in terms of a normalized pure state \tp(t)} which follows a stochastic "trajectory" in 
Hilbert space. By averaging the pure state projector \ip}(tp\ over all possible trajectories with 
appropriate weights, one reproduces the density operator p = M (\ip) (ip\) . This is analogous 
classically to replacing the Fokker-Planck equation for probability densities with a stochastic 
Langevin equation for single trajectories. 

Unfortunately, unlike the case of classical Brownian motion, the unraveling of the master 
equation (|l|) is not unique. Thus, there is some ambiguity in how one separates classical and 
quantum uncertainties, related to the ambiguity in identifying density matrices with quantum 
ensembles. In this section we consider two well known unravelings of (|l|). 

In quantum state diffusion (QSD), the (Ito) stochastic evolution equation for the normalized 
state vector \ip(t)) reads: 



|#(*)> = -iH\m)dt -\Y.( L ) L i ~ 2 ( L ]h L j + \(L,)^\ 2 )\m)dt 

j 

+ £(L i -(L i )vOIV>(*Mj (2) 



j 



where the "noises" d£j are complex- valued Wiener processes of zero mean M(d^j) = and 
correlations M{d^jd^k) = 0, M(d^d^k) = Sjkdt. This equation describes a continuous non- 
differentiable evolution similar to the familiar diffusive paths of a classical Brownian particle, 
but in Hilbert space instead of real space. QSD is the only continuous unraveling which satisfies 
the same symmetries as the master equation itself jy], |3| . 

Our second example is the quantum jumps (QJ) unraveling, which is closely related to 
photon counting [12, 13]. However, it can be defined for any master equation (|l|). The stochastic 
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increment for the wave-function is 



j 

• L 3 \m) 

ssume the 

the evolution is continuous and differentiable. However, whenever dNj = 1 there is a "jump 



+i: ,j^m_ m)]iNj (3) 



The discrete Poissonian noises dNj assume the values or 1. Most of the time dNj = and 



to the state Lj\ip(t)} / J (L^Lj) t - The dNj processes have mean values M^(dNj) = (L^L^^dt 
and correlations dNjdt = and dNjdN^ = Sj^dNj. This means, essentially, that jumps occur 
randomly with an average rate (L^Lj). 

It has been shown that a quantum jump description of balanced heterodyne detection pro- 
duces an equation identical to the QSD equation (^) in the limit of a strong local oscillator 



[14, 15], and that a QSD description of photon detection produces jump-like behavior [O]. 
Also, in the classical limit (i.e., high photon numbers), the quantum jump equation begins to 
exhibit diffusive behavior as well |16|| . Thus, though these unravelings are distinct, they can 
behave similarly in certain limits. 



3 Harmonic oscillator 

Let us illustrate these two unravelings for a simple example: the damped harmonic oscillator 
at finite temperature. H = uia^a, L\ = ^/n/y and L2 = y/(n + 1)7 a where n is the thermal 
equilibrium mean photon number, n = {a)a) p , and 7 is the inverse relaxation time. For QSD one 
can show that any initial states tends to a coherent state: \ip{t)) — * \ctt), where a\at) = ott\at) 
and at = ((q)ip+i(p) 1 / J )/V2. Furthermore, the evolution of at is governed by a classical equation, 

dat = —iuatdt — —atdt + \/nrjd£t (4) 

Hence, for this example at least, the QSD equation fully describes how the environment localizes 
the quantum state down to a minimum Gaussian wavepacket, and how this wavepacket follows 



a classical trajectory |17|, |18| |19|, |20| |2l|, |22|. 

For quantum jumps a similar analysis can be carried out, and for a system in a coherent 
state the evolution of at is governed by the equation 

7 (X 

dat = —iujatdt atdt + yrry^p dWt , (5) 

2 V M + 1 

where dWt is a real, non-Gaussian stochastic differential variable with mean M{dWt) = and 
M{dWtdWt) = dt, representing the difference between the mean of the jump processes and their 
actual values. In the classical limit of large a, this dWt approaches a diffusion process, and 
equation (||) has behavior essentially identical to (Q). 

Somewhat surprisingly, states tend towards small wavepackets for quantum jumps as well as 
QSD, so the same "classical limit" seems to exist for both [J16|j . 
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4 Chaotic Duffing oscillator 



A good example of a nonlinear system is the forced, damped Duffing oscillator. This has a 
classical equation of motion 

— + 2T— + x 3 -x = gcos(t), (6) 

and for some choices of T and g is chaotic [^] . 

Because the equation of motion includes explicit time-dependence, the solutions lie in a 
three-dimensional phase space x,p,t. It is helpful to consider a discrete surface of section of 
this system. Let (xo,po) be the initial point of the forced, damped Duffing oscillator at time 
to = 0. Then we can define a constant phase map in the x-p plane by the sequence of points 
(x n ,p n ) = (x(t n ),p(t n )) at times t n = 2irn. Figure 1 illustrates this in the chaotic regime, where 
we can clearly see from the surface of section the fractal structure of the strange attractor. 

Quantizing the Duffing oscillator is straightforward using the QSD formalism. The Hamil- 
tonian operator is 

H(Q,P,t) = P 2 /2m + Q 4 /4-Q 2 /2 + gcos(t)Q + VT{QP + PQ), (7) 

and the damping is represented by a Lindblad operator 

L = 2Vfa = V2T(Q + iP), (8) 

where we have assumed h = 1. The last term in the Hamiltonian is an ansatz, added to give 
the correct equations of motion in the classical limit; it is necessary due to the simplistic model 
of the dissipative environment. 

This system is far from classical. To go to the classical limit, we introduce a scaling factor 

P, 

Hp(Q, P, t) = P 2 /2m + Q 4 /4f3 2 - Q 2 /2 + g(3cos{t)Q + Vf(QP + PQ). (9) 

As we increase /3, the scale of the problem (compared to Ti) increases by (3 in x and p, without 
altering the dynamics. Thus, (3 —* oo is the classical limit of this system. Classical behavior 
should emerge from the system in this limit. This is supported by the numerical calculations. 
(See figures 1 and 2.) 

The classical unsealed problem is bounded within a small region of phase space. (See figure 
1.) Clearly, the quantized version of this problem (with h = 1) should be far from the classical 
limit. One would expect to observe little trace of the classical fractal structure. The expectation 
values of Q and P should be dominated by noise. We solved the forced, damped Duffing oscillator 
numerically by integrating the QSD equation for three different values of the parameter /3. 

Examining figure 2a, we see that the results of the numerical calculation match our expec- 
tations very well. The expectation values appear randomly distributed; they are dominated 
almost completely by the stochastic terms of the equation. 

Approaching the classical limit of large /3, more and more of the classical structure of the 
attractor appears. At first the broad outlines of the attractor are formed, then increasing levels of 
substructure. (Figures 2b-2c.) The full fractal structure of a strange attractor is only attainable 
in an unattainable purely classical limit. In an actual physical system, the uncertainty principle 
provides a lower cutoff to the scale- invariance of the strange attractor. 

It can easily be shown that the quantum jumps formalism has a "diffusive limit" similar to 



QSD as one approaches the classical limit [16]. Thus, exactly the same qualitative behavior is 
demonstrated in that case as well. 
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5 Conclusions 



Neither QSD nor QJ can be regarded as a fundamental theory, since both divide the universe 
arbitrarily into a system and environment; and they represent two different and interesting un- 
ravelings of the master equation (though only QSD has the same symmetries as the master 
equation [pd| , This is similar to the ambiguity which arises in consistent histories, in which 
it is possible to choose different sets of histories corresponding to different coarse-grained de- 
scriptions 24]. On the other hand, QSD (or QJ) makes testable predictions, and demonstrates 
the rise of classical physics in a way which is difficult or impossible to see with other theories. 

Without entering into a sterile debate, let us note that there seems to be an interesting link 
between consistent histories and the stochastic models used in this contribution. It appears that 
the solutions of the stochastic equations provide consistent sets of histories. Hence, it seems 
that the two approaches are related in a way similar to the relationship in classical mechanics 
between the global view given by the Maupertuis principle and the local view given by the 
Newton equations. The global view may be more elegant, but the local one is more useful to 
compute actual numbers J2f| p6| . 

Quantum mechanics is nonlocal. Classical mechanics is local. Consequently classical me- 
chanics can not explain all quantum phenomena. Conversely, it is cumbersome to use quantum 
mechanics to describe classical phenomena. Not only are the computations more complex, but — 
and this is the main point — it is conceptually more difficult: one has to argue that nonlocality, 
entanglement and the principle of superposition can be set aside when crossing the "quantum — ► 
classical" border. Clearly, nonlocality, entanglement and the principle of superposition should 
become irrelevant in the classical limit. But why should one argue? Shouldn't it just come 
out of the equations? Does it come out of the equations? This contribution is about the last 
question. And the answer is: "it depends on which equation." 
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Figure 1. The constant phase surface of section for the classical forced, damped Duffing 
oscillator in the chaotic regime, T = 0.125, g = 0.3. 

Figure 2. The constant phase surface of section for a single QSD trajectory of the quantum 
forced, damped Duffing oscillator in the chaotic regime, T = 0.125, g = 0.3, for three scalings: 
a) p = 1, b) = 4 and c) p = 10. 
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Forced, Damped Duffing Oscillator, Chaotic Regime (QSD) 
Surface of Section, beta = 1, 4, 10 
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